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This talk provides a summary of the resuhs obtained by the HotQCD collaboration on the equation of state and 
the crossover transition in 2+1 flavor QCD. We investigate built thermodynamic quantities - energy density, 
pressure, entropy density, and the speed of sound over the temperature range 140 < T < 540 MeV. These 
results have been obtained on lattices of temporal size Nr = 6 and 8 and with two improved staggered fermion 
actions, asqtad and p4. Our most extensive results are with masses of the two degenerate light quarks set 
at mi = 0.1ms corresponding to the Goldstone pion mass m-,^ between 220 — 260MeV. In these simulations, 
the strange quark mass is tuned to its physical value and constant values of mi/ms define lines of constant 
physics. We also summarize the current state of results on observables sensitive to the chiral and deconfining 
physics — the light and strange quark number susceptibilities, the chiral condensate and its susceptibility, and 
the renormalized Polyalcov loop. Our results indicate that the deconfinement and chiral symmetry restoration 
occur in the same narrow temperature interval. 



1. Introduction 

Ongoing experiments at RHIC and proposed exper- 
iments at LHC aim to understand the properties of 
hot dense nuclear matter created in the collision of 
two relativisitic nuclei. At sufficiently high tempera- 
tures and densities RHIC data support the creation 
of a quark-gluon plasma in the central region that 
undergoes a transition back to hadronic matter as it 
expands and cools. The goal is to explain the creation 
and evoluton of this medium. Hydrodynamic descrip- 
tions used to model this evolution provide a good fit 
to the data and are thus the phenomcnological tool of 
choice. First principle calculations using lattice QCD 
yield a number of properties of QCD as a function 
of temperature that are essential inputs in these hy- 
drodynamical analyses. These properties include the 
nature of the transition (with respect to both con- 
finement and chiral symmetry breaking) between the 
quark-gluon plasma and hadronic matter, the transi- 
tion temperature, the equation of state of QCD and 
transport coefficients as a function of temperature in 
the range 140 - 700 MeV. 

HotQCD is a US wide collaboration engaged in the 
study of QCD at finite temperature and density using 
lattice QCD ^. It brought together members of the 
MILC and RBC-Bielefeld collaborations to carry out 
large scale simulations on IBM Bluegene/L supercom- 
puters at the Lawrence Livermore National Lab and 
on the NYBlue at the New York Center for Compu- 
tational Sciences at BNL. 

Our goals are to perform detailed simulations of 
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2-1-1 flavor QCD using staggered fermions at the physi- 
cal values of the strange and light quark masses. Most 
of the results presented here are for mz = O.lm^, about 
a factor of two heavier than the mean physical u and 
d quark mass, mi O.OSSm,,. Improvements in algo- 
rithms and computer resources are being used to in- 
crementally approach the physical up and down quark 
masses and the continuum limit, thus providing high 
precision results. Also, in the current simulations the 
up and down quark masses are taken to be degenerate. 
The physical quantities we are calculating include: 

• The Equation of State (EoS) of QCD over the 
temperature range 140 — 700 MeV that is being 
probed by relativistic heavy ion experiments at 
Brookhavcn and will be studied in more detail 
in the future at the LHC. 

• The nature of the deconfinement transition be- 
tween the hadronic phase at low temperature 
and the quark-gluon plasma at high tempera- 
tures. All simulations with staggered fermions 
show a rapid crossover rather than a genuine 
phase transition for physical values of the light 
quark masses. In the absence of a phase transi- 
tion, there is a priori no unique transition tem- 
perature as it can depend on the probe. Thus, 
the temperature at which this transition takes 
place remains a subject of investigations. The 
status of current estimates is discussed at the 
end of this paper. 

• The restoration of chiral symmetry at high tem- 
perature and whether this chiral transition is co- 
incident with the deconfining transition. 

• A detailed understanding of the physics in the 
transition region and the approach of thermo- 
dynamics quantities such as pressure, entropy, 
energy density and the speed of sound, to the 
Stefan-Boltzmann limit. 
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To understand and control systematic errors in the 
lattice calculations we are taking the following steps: 

• Simulations at Nr = 4, 6, and 8 have been per- 
formed with two improved staggered fermion ac- 
tions - asqtad and p4 - to understand O(a^) 
discretization errors. 

• Extrapolation to the continuum limit to remove 
discretization errors will be made using simula- 
tions at Nr = 6,8 and a new set of = 12 
lattices. 

• The theory possesses a chiral phase transition in 
the limit to; = 0. To study the chiral behavior, 
extrapolations to the physical light quark mass 
limit and the chiral limit will be made using sim- 
ulations at mi/rris = 0.2, 0.1 and 0.05. 

• The transition region, 160 — 230 MeV, is being 
sampled in intervals of 5 MeV or less using high 
statistics simulations in order to pin down the 
position of the peak in the chiral susceptibility 
and the inflection point in the light quark num- 
ber susceptibility. 

The first part of this report summarizes results for 
the equation of state obtained by the HotQCD col- 
laboration In the second part we discuss our re- 
sults for the temperature dependence of the chiral and 
deconfining transitions and finally compare these re- 
sults against recent data presented by the Wuppertal- 
Budapest collaboration 0. 

2. Setup of Lattice Calculations 

To determine the lattice parameters for a systematic 
finite temperature analysis, zero temperature simula- 
tions are used to first set the lattice scale a and the 
physical value of the strange quark mass as a func- 
tion of p. The scale a is extracted using the Sommer 
scale To obtained from the derivative of the static qq 
potential evaluated on the lattice 

\ / r—ro 

with ro = 0.469(7) fm taken from the Upsilon 2S-IS 
splitting [3]. The strange quark mass is then set using 
the ss pseudoscalar meson state with Mssro — 1.58. 

The last quantity to fix are the light (up and down) 
quark masses that are taken to be degenerate. We de- 
fine lines of constant physics (LCP) by keeping mi/mg 
fixed. In our calculations three values of mi/ms ~ 
0.2, 0.1 and 0.05 are being simulated. Results for the 
physical value of the light quark mass are then ob- 
tained by extrapolation in to; at fixed lattice scale. 
With these definitions, the continuum limit along a 



LCP can be taken by varying a single paramter, the 
gauge coupling /3 ^ oo or equivalently a since 
P and a are related by dimensional transmutation in 
QCD. The results presented here are obtained mostly 
for mi/nis = 0.1 with Nr = 6 and 8 lattices, so es- 
timates of extrapolated values and uncertainties in 
them should be considered preliminary. 

In the finite temperature calculations, the tempera- 
ture T, lattice scale a and the size of the lattice in the 
temporal direction Nr are related as T = 1 /aNr . For 
fixed Nr we increase T by decreasing a or equivalently 
increasing /3, and adjusting mg for each /3. 

There are two overarching issues with the use of 
staggered fermions. Both have to do with the fact that 
staggered fermions preserve a residual chiral symme- 
try at the expense of a 4-fold doubling of flavors, i.e., 
every continuum quark flavor develops four copies on 
the lattice. These four copies are called taste. The 
first issue with using staggered quarks is one of prin- 
cipal. To simulate a single flavor (say an s quark) 
we take the fourth root of the fermion determinant 
in calculating the fermion force in Monte Carlo sim- 
ulations. Creutz ^ has argued that a strong mix- 
ing of tastes with different chiralities leads to an in- 
correct 't Hooft vertex, and concludes that rooting 
can often be a good approximation but predictions 
for non-perturbative physics where the 't Hooft vertex 
is important can not be trusted. On the other hand 
Sharpe [5] and Goltermann ^ argue that while this 
rooting trick is ugly a well defined continuum limit 
exists and the theory recovered in this limit is QCD. 
Numerical results obtained with staggered fermions 
when compared to experimental data indicate that the 
rooting trick works i.e. this circumstantial evidence 
provide support that staggered simulations reproduce 
QCD. Clearly, our calculations rely on the validity of 
the arguments reviewed by Sharpe and Golterman. In 
short, our work does not shed any independent light 
on this issue. 

The second issue is that, at finite lattice spacings, 
the 4 copies called tastes, are not degenerate and the 
taste symmetry is badly broken for both asqtad and p4 
actions at 1/a < 2 GeV corresponding to the transi- 
tion region for Nr < 8 simulations. This issue is being 
explored and improved versions of the staggered ac- 
tion are being designed to reduce this breaking. The 
most severe impact of this breaking is in the pion sec- 
tor. In our p4 action calculations with mi/rris = 0.1 
on Nr = 8 lattices, the mass of the Goldstone pion is 
M^ro ~ 0.52 corresponding to « 220 MeV while 
the other fifteen higher mass taste flavor pions are 
at 440 MeV or heavier. The splitting for asqtad ac- 
tion is slightly smaller. We expect this violation of 
taste symmetry to affect the results in the low tem- 
perature confined phase and in the transition region 
where the dominant excitations are pions. Unfortu- 
nately, at this point we do not have a good estimate 
of the uncertainty introduced in the various quantities 
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we calculate due to taste symmetry breaking. 

The single basic quantity given by lattice simula- 
tions en-route to calculating thermodynamics vari- 
ables such as the pressure and energy density is the 
trace anomaly Q^^'^a* — {e — 3p)a^. We first discuss 
the data for it and the extraction of the EoS from 
it. The second part introduces probes used to under- 
stand the nature and location of the deconfining and 
chiral symmetry restoration transitions and presents 
our data. Finally we compare these results to those 
from the Wuppertal-Budapest collaboration [5]. 



3. Trace Anomaly (e - 3p)/T^ 

The calculation of G^^a* involves calculating terms 
that make up the gauge and fermion actions. As a re- 
sult, most of the computer time is spent in generating 
the ensemble of statistically independent gauge con- 
figurations. Results presented here represent over 100 
million node hours on the BlueGene/L. 

Data for (e — 3p)a^ is converted into physical units 
using the scale tq and subsequently in units of the 
temperature T as shown in Fig. [l] Overall, the data 
show that the asqtad and p4 actions give consistent 
results, with differences consistent with the expected 
magnitude of O(a^) errors. The details are highlighted 
in figures [2] [3] [4] and discussed below. 

From the trace anomaly the pressure is obtained 
using the relations 
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The energy density, entropy density {e+p)/T'^ and the 
speed of sound are then given by appropriate combi- 
nations of (e — 3p)/r'*and p/T'^. 
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Figure 1: Data for the trace anomaly (e — 3p)/r* calcu- 
lated using the asqtad and p4 actions at A'',- = 6, 8. This is 
the single basic quantity obtained from lattice simulations. 
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Figure 2: Details of (e - 'ip)/T^ in the range 140 - 200 
MeV. 



Low Temperature Region: Comparing Nr = 6 
and 8 data in Fig. [2] for the low temperature region 
show a ~ 5 MeV shift towards lower temperatures 
under Nt = 6^8. This shift is indicative of 0{a?) 
discretization errors. To calculate the pressure using 
Eq. |2] we need to fit the data obtained at a discrete 
set of points to a functional form that can be inte- 
grated and fix the normalization of at some low 
starting temperature Tq. We first discuss the normal- 
ization and return to the determination of the func- 
tional form in the next section. For T < 150 MeV, 
the discretization errors in lattice simulations, espe- 
cially those due to the large taste symmetry violations 
in asqtad and p4 staggered actions, are pronounced 
so Nr = 8 data have unknown systematic uncertain- 
ties. We, therefore, explore predictions of the hadron 
resonance gas (HRG) model as it is an independent 
method in which control over systematic uncertain- 
ties improves as T is taken below the pion mass scale 
= 140 MeV [g. In Fig. [2] we plot results for 
two choices of the cutoff scale rrimax up to which res- 
onances are included, nimax = 1.5 GeV (dot dashed 
black curve) and 2.5 GeV (upper dashed curve). Un- 
fortunately, the HRG estimates show significant sensi- 
tivity to the cutoff rrijnax for T > 100 MeV. Also, there 
is no matching point Tq between lattice and the HRG 
data. There is a trend in the lattice results to move 
towards the HRG values as a —> 0, i.e., as Nr = 6 — > 8, 
however, at his point, we do not have a way to quan- 
tify it, and conclude that further work is required to 
reduce the uncertainty associated with the choice of 
normalization. 

Intermediate (Peak) Region: Figure |3] shows 
that the peak in e'^'^/T'* occurs at T > 200 MeV. 
The differences between the asqtad and p4 results for 
Nr = 6 versus 8 are the largest in the interval 200 < 
T < 300 and show a 15 - 20% effect at the peak. We 
also find that the peak in the asqtad data is shallower 
and shows smaller difference between Nr = 6 and 8 
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values as compared to the p4 results. 
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Figure 3: Details of (e - 3p)/T^ over the range 180 - 300 
MeV. 

High temperature region T > 300 MeV: Data 
in Fig. |4] show that the two actions give consistent 
results and the difference between Nr = 6 and 8 are 
small. Q^^^ /T^ should, at sufficiently high tempera- 
tures, vanish as g^{T), however, between 300 — 700 
MeV the data show a much faster variation. We, 
therefore, use the ansatz 
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(3) 



to fit the data. Fits show no sensitivity to the leading 
term. Consequently, this term is neglected in the 
fits shown in Fig. [4] and used for the extraction of p. 

To improve these fits and determine the shape of 
the curves for T > 300 MeV two enhancements to our 
data sets are needed. First, we need data at more 
points. Second, checks for finite volume corrections 
need to be made for T > 500 MeV when the lat- 
tice size in the spatial directions used in the simu- 
lations, Ns — 32, approaches Nga — 1/160 for a in 
MeV~^. Here a corresponds to the zero temperature 
scale at each of the the gauge couplings used to simu- 
late T > 500 MeV lattices, and the temperature scale 
160 MeV is used as it marks the rapid onset of thermal 
fluctuations as shown in Fig[l] 

3.1 . Fits to and extraction of p 

To extract the behavior oip versus T using Eq. [2] we 
need to fit the data for Q'^^/T* calculated at a discrete 
set of points with a smooth function. We do this using 
a piecewise smooth interpolation. For T < 170 MeV 
an exponential fit is used, while for 170 < T < 250 
MeV separate quadratic fits over several intervals were 
used. We found it necessary to match the function 
and its derivative at each end of these intervals to get 
a smooth curve for p. For T > 250 MeV the result of 
the fit with the form given in Eq. |3]is used. These fits 
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Figure 4: Details of (e - 3p)/T^ over the range 300 - 700 
MeV. 



are shown for both actions and for = 6 and 8 in 
Figures |5] 

The fits are very similar in shape for the two ac- 
tions. In each case shifting the fits to Nr = 6 data 
by « 5 MeV to lower temperatures give much better 
agreement except for the height of the peak in the p4 
data. This shift by 5 MeV suggests the size of dis- 
cretization errors coming from a combination of (i) 
the errors in the finite T simulations and (ii) from the 
scale a extracted from zero temperature simulations 
clustered around (3 values corresponding to the two 
peak regions for Nr — 6 and 8. 

Our results for p using Eq. [2] (and consequently all 
other thermodynamic quantities) are obtained with 
the normalization p = OaiT=TQ = and requiring 
the fit to the lattice data to pass through this point. 
Redoing the fits starting with the HRG value at T = 
100 MeV results in a global shift in the pressure and 
energy density curves by sa 0.8. The size of this effect 
is illustrated by the black square at T = 550 MeV in 

Once (e — 3p)/T^ and p/T^ are known we get the 
other thermodynamics quantities; entropy density as 
(e + p) /T^ and the speed of sound by using the ratio 
p/e in 



(4) 



2 _ dp _ d{p/e) p 
" de " ' de ^ e ■ 

These results are shown in Figures [7] and [Sj We regard 
the spread in the results between the two actions as 
indicative of discretization errors. 

In Figures [6] and [7] we highlight the band 185 — 195 
MeV which we estimate covers the inflection point in 
the energy density e and the entropy density data. 
Both of these thermodynamics quantities probe de- 
confinement — the change in the thermodynamics as 
the degrees of freedom evolve from hadrons to quarks 
and gluons due to the change in their free energies. 
Note that while we discuss other probes of decon- 
fincmcnt and chiral symmetrey restoration later, it is 
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Figure 6: Results for the energy density e/T'* and the 
pressure 3p/T^ for the two actions obtained on — 8 
lattice with m/ms = 0.1. 



these thermodynamic quantities that are inputs in the 
hydrodynamic analysis of RHIC data. 

Our data show that for T > 300 MeV the degrees 
of freedom are quarks and gluons and the deviations 
from the Stefan-Boltzmann Hmit are 10 — 15%. Thus 
experiments at the LHC, which will probe the full 
range 140 < T < 700 MeV, will be able to explore 
the behavior and nature of collective excitations in 
the quark-gluon plasma. Even though the HotQCD 
estimates for the EoS do not include extrapolation to 
the continuum limit and physical light quark masses, 
we propose that they are realistic and ready for use in 
phenomenological analyses. 



4. The Deconfining and Chiral Transition 

The data presented in previous sections show that 
the transition from hadronic matter to quark-gluon 
plasma is a rapid crossover and not a genuine phase 
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Figure 7: The entropy density obtained on lattices with 
temporal extent A'^^ = 6 [51 [TO] and 8. 

transition. This is consistent with all existing calcu- 
lations with staggered fermions. We expect a singu- 
larity in the partition function to exist in the chiral 
limit. What remains to be settled for QCD with a 
physical value of the strange quark mass is whether 
this phase transition occurs only at zero light quark 
masses (m; = merit = 0) or at small but nonzero 
value (m; = merit > 0). In the first case the tran- 
sition would belong to the universality class of the 
3-dimensional 0(4) symmetric spin model and in the 
latter case to the 3 — d Ising model jTT] . 

In light of this expected behavior, the questions that 
we would like to answer regarding the nature and loca- 
tion of the transition from hadronic matter to quark- 
gluon plasma in QCD are: 

• How is the observed crossover in thermodynamic 
quantities for physical light and strange quark 
masses influenced by the chiral phase transition? 

• What do standard probes of deconfinement (ex- 
pectation value of Polyakov loop and quark 
number susceptibilities) and chiral symmetry 
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Figure 8: Data for p/e and the speed of sound extracted 
from it. Lines without data points are obtained using the 
interpolating functions shown in Figures |5] The dash-dot 
hne is the HRG estimate. 



Figure 9: Data for the Polyakov Loop {Lren}- Note the 
gradual rise in contrast to the rapid crossover seen in ther- 
modynamic quantities. 



restoration (chiral condensate and the chiral sus- 
ceptibility) tell us? 

• Do the crossovers indicating deconfinement and 
chiral symmetry resoration happen at the same 
temperature (we will denote these as Tc but note 
that it represents the crossover temperature and 
not a critical temperature) and what are the 
mechanisms that drive them? 

Our results for the various probes are presented next. 



4.1. Renormalized Polyakov Loop 

The logarithm of the expectation value of the 
Polyakov loop is related to the free energy of an iso- 
lated quark, 



{L)^exp{~FjT) 



(5) 



In QCD (L) is not a genuine order parameter, i.e., zero 
in the confined phase and non-zero in the deconfined 
phase, for any finite values of the quark masses. As 
far as we know, the Polyakov loop is not directly sen- 
sitive to the singular structure of the partition func- 
tion in the chiral limit. It does , however, exhibit a 
broad crossover as shown in Fig. |9] where we plot the 
renormahzed quantity (Lren) = Z{[3)^^ (L) [TU]. The 
renormalization is needed to remove the self-energy 
contributions to the free energy. 

Data for (Lren) exhibit a gradual rise rather than 
a rapid crossover. This feature is also reflected in the 
lack of a peak in its susceptibility. We find only a 
broad shoulder after the rise for all data with N^. > 6 
lattices. Since there is no well defined infiection point 
in (Lren) Or a peak in its susceptibility we do not 
consider it a reliable probe and argue it should not be 
used in precision studies to determine the location of 
a crossover Tr indicative of deconfinement. 



4.2. Quark Number Susceptibility 

Thermal fiuctuations of the degrees of freedom that 
carry a net number of light or strange quarks are 
probed by taking derivatives of the partition function 
with respect to the corresponding chemical potential. 
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=0 



The quark number operator ipjo'^P is a good probe in 
lattice calculations since it does not require normaliza- 



tion. It counts the charge 



the net number of 



(N^). There- 



light or strange quarks. Similarly Xq 
fore, in the continuum and infinite temperature limit 
limT^oo Xq/T'^ = the Stefan-Boltzmann value for 
an ideal massless one flavor quark gas. The approach 
to unity in our data for the light quarks is shown in 
Fig. [To] 
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Figure 10: Data for the light quark number susceptibility 
Xi/T^- Note the crossover is consistent with the band at 
185 — 195 MeV and the value approaches unity rapidly 
after 250 MeV. 
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As the transition region is approached from above, 
hadronic states with appropriate quantum numbers 
start to contribute in proportion to their Boltzmann 
weight ~ exp(— A//T). At low temperatures the dom- 
inant contributions are from the hghtest hadronic 
state that has the correct quantum numbers, i.e., 
Xs/T^ - eM-Mx/T), while Xi/T\-^ exp{-MjT). 
Since Mk remains finite in the limit m; — > 0, only 
Xi /T^ is directly sensitive to singularities in the QCD 
partition function in the chiral limit. This expecta- 
tion is confirmed by c omp aring the behavior of e/x^T^ 
versus e/xiT'^ in Fig- H Below the transition region 



e/XsT'^ rises rapidly in contrast to e/xiT'^- In fact it 
will diverge because Xs is not sensitive to the lightest 
hadrons and goes to zero faster than e which does re- 
ceive contributions from the pions. We therefore claim 
that xi/T^ and not Xs/T^ should be used to probe the 
singular structure of the theory. 



Both renormalizations are performed by defining the 
following subtracted quantity for light quarks 



Ai,.(T) 



LT 



(V''0)/,O 



(8) 



The subtractions proportional to (■0V') s remove the 
singular part of the additive renormalization propor- 
tional to niq/a^. Division by the T = estimate ob- 
tained at the same lattice scale a cancels the multi- 
plictive renormalization as it is T independent. 

Data for A/^<;(T) along the LCP given by m; = 
0.1ms are shown in Fig. [T2j We find a rapid crossover 
in the same temperature range as exhibited by the 
bulk thermodynamic observables and the light quark 
number susceptibility. This indicates that deconfine- 
ment and chiral symmetry restoration happen at the 
same temperature. 
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Figure 11: Comparison of the strange and light quark 
number susceptibility in units of the energy density, i.e. 
e/XsT^ versus e/xiT^- e/Xs^^ is seen to grow rapidly and 
is expected to diverge as mi ^ since e gets contributions 
from the pions while Xs as exp(— A/k/T). 

The data for Xi/^^ in the transition region (see 
Fig. 10 1 show a rapid crossover with the location of 



the infiection point again consistent with the band 
185 — 195 MeV. We therefore conclude that all probes 
of the deconfinement transition, £, the entropy density 
and xi/T"^: give consistent results. 

4.3. Chiral Condensate and 
Susceptibility 

The light and strange quark chiral condensates are 
given by the derivatives of the free energy with respect 
to the quark masses: 



T dlnZ 
V dnia 



q^l, 



(7) 



Figure 12: p4 and asqtad data for the renormalized chiral 
condensate Ai^s{T) along the LCP given by mi = 0.1ms. 
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Figure 13: Same data as in Fig. [12] for the renormalized 
chiral condensate A;,s(r) but with the A*",- = 6 data shifted 
to the left by 5 MeV for each action. 



At finite values of the quark masses {'>pip)q requires 
both multiplicative and additive renormalizations. 



To estimate the magnitude of the difference between 
Nr ~ 6 and 8 data we plot, in Fig. [13] the same data 
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after shifting the Nt- = 6 values by 5 MeV to the 
left. The results for = 6 and 8 lattices now come 
together. The replot also brings out the difference in 

results between the two actions the difference is 

in the shape of curves. The data with the p4 action 



shows a slightly steeper crossover. Based on Fig. 13 
we estimate that both systematic effects are ~ 5 MeV 
and the location of the inflection point is covered by 
the band 185 — 195 MeV shown by the vertical lines. 

The chiral susceptibility is given by the second 
derivative of In Z with respect to the quark mass m: 



Xl,s 

rp2 



T ( d'^lnZ 



V \d^mi.s/T)^ 
(Tr(AC)) 

+ I (((TrM,;;i)2 



(TrAf, 



-1\2 



(9) 



where the first term in the right hand side gives the 
contributions of connected Feynman diagrams and the 
second gives those from the disconnected diagrams, xi 
is the best probe of the chiral transition as it exhibits 
a peak even when the theory has a crossover and it 
will diverge at the chiral phase transition, i.e., in the 
chiral limit. The location of the peak is used to mark 

Preliminary unpublished data from the N-r = 8 sim- 
ulations is shown in Fig. [14] The left hand figure 
shows data for mi/nis = 0.2, 0.1 and 0.05 for the two 
actions. We find that the height of the peak grows 
with decreasing quark mass and is expected to diverge 
with an exponent 1/^—1 corresponding to the 0(4) 
unversality class of the chiral transition. 

The second noteworthy feature is that the peak is 
broad and its width increases with decreasing quark 
mass. This is also expected because in the chiral limit 
the pions are massless for all values of T < Tc and 
give rise to a singularity in the "magnetic" part of 
the free energy, and the chiral susceptibility is ex- 

1 /2 

pected to diverge as ' [12^ . The plot on the right 
hand side shows the same data scaled by (mi/ms)^^^ . 
The peaks in the three data sets collapse together, 

1/2 

consistent with the expected m/ behavior. We do 
not find simultaneous evidence for the expected 0{N) 
scaling of the height of the peak mentioned above. 
This dual scaling behavior raises the question — what 
point should be labeled as Tc? In the absence of a 
clear peak, we will choose the right edge to denote 
the chiral symmetry restoration temperature Tf^"^^. 

The third feature shown by the data is that the 
location of the peak moves to lower temperatures with 
decreasing quark mass. For m; = 0.1ms the location is 
consistent with the 185—195 MeV range found earlier, 
however, the trend suggests a lower value (by 5 — 10 
MeV) in the physical quark mass limit [T^] . Similarly, 
the location of the peak in the asqtad action is shifted 
with respect to p4 to lower temperatures at each value 
of the quark mass by ~ 5 MeV. 



The bottom line is that we have identified a number 
of systematic uncertainties that are each ~ 5 MeV and 
go in the same direction of lowering the estimate for 
Tc. It is therefore necessary that the continuum along 
nil ~ 0.038 LCP is taken in order to extract physical 
results. 



5. Continuum Limit and HotQCD Plans 
for the Future 

The results presented in the previous sections have 
been obtained on lattices with Nr = 6 and 8 along a 
LCP defined by mi — 0.1ms. Two extrapolations are 
needed to obtain physical results. These are extrap- 
olations in the light quark mass to the physical value 
mi w 0.038 and the continuum limit — s- oo with 
a ~f 0. The present data do not allow us to make 
either extrapolation reliably. 

A number of simulations are in progress to address 
these two limitations: 

• For both asqtad and p4 actions and for Nr = 
6 and 8 lattices we are completing simulations 
at three values of the light quark masses, i.e. 
mi/mg = 0.2, 0.1 and 0.05. These data will 
allow us to make extrapolations to the physical 
quark mass and the chiral limit. 

• Having established in this study that asqtad 
and p4 actions give consistent results, we have 
started new asqtad simulations with = 12 
and mi = 0.05ms. The combined data at Nr = 
6, 8 and 12 will allow us to take the continuum 
limit along the LCP defined by mi = 0.05ms 
and extrapolated data at m; = 0.038ms. 

• We are exploring the Highly Improved Staggered 
Quark (HISQ) action [T?] which has been de- 
signed to improve both the O(a^) scaling behav- 
ior and reduce taste symmetry violations [i TS] . 
As mentioned earlier, the uncertainties associ- 
ated with taste symmetry violations are not well 
understood and could be large. 



6. Comparison of with Results from 
the Wuppertal-Budapest Collaboration 

The Wuppertal-Budapest collaboration (WBC) has 
recently presented the following results on the location 
of the crossover indicating deconfinement and chiral 
symmetry restoration: 



Tc = 146(2) (3) MeV Chiral, 
Tc = 170(4) MeV Deconf. 



(10) 



WBC results include extrapolation to the continuum 
limit using Nr = 8, 10,12 lattices along the m/ = 
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Figure 14: Preliminary HotQCD data for the disconnected part of the chiral susceptibility Xi{T) for the p4 and asqtad 
actions with Nr = 8 and three values of the quark mass mi/ms — 0.2, 0.1 and 0.05. The right panel shows the same 
data as in the left panel but scaled by [mi/msY^^ ■ The behavior of the peak as a function of the quark mass is discussed 
in the text for both panels. 



0.04tos LCP and are therefore labelled "physical". 
They raise the following two issues when comparing 
to the HotQCD results discussed earlier. 

• The WBC estimate for the chiral transition at 
146(2) (3) MeV is about 40 MeV lower than the 
HotQCD results. 

• The crossover in the chiral symmetry and de- 
confinement sectors are distinct and separated 
by « 24 MeV whereas all HotQCD data indi- 
cate they are coincident. 

Of these, we consider the difference in the location 
of the chiral transition the more significant issue that 
needs to be resolved. The reasons for discounting the 
difference in the location of T^'^^°^^ as significant are 

• The best HotQCD data are on Nt ~ 8 with 
mi = 0.1ms lattices. Three systematic effects, 
each shifting the estimate by « 5 MeV to lower 
temperatures, have been identified. These are 
(i) the difference between Nr — 6 and 8 lattice 
data; (ii) extrapolation in light quark masses 
and (ill) the difference between asqtad and p4 
data. A fourth possible shift of 5 MeV arises 
from the recent change in the experimental value 
of fii which resulted in lowering the WBC re- 
sults by ~ 5 MeV [5]. This change affects 
the continuum extrapolation carried out by the 
WBC and not the HotQCD results presented at 
fixed a with the scale set by rg. All four fac- 
tors are not necessarily independent, however, 
the point we want to make is that the difference 



between 170(4) and 185 - 195 MeV can largely 
be accounted for by a combination of these sys- 
tematic effects. 

• More important, we have presented reasons for 
why the two methods used by the WBC to ex- 
tract T^econf^ ^YiQ Polyakov loop and strange 
quark number susceptibility, do not probe the 
singular structure of the theory. Our contention 
is that all previous results based on these probes 
do not qualify as precision measurements. 

• Probes sensitive to the singular structure (en- 
ergy and entropy density and the light quark 
number susceptibility) show the two transitions 
to be coincident for all values of the lattice pa- 
rameters simulated by the HotQCD collabora- 
tion. It remains to be seen if this conclusion sur- 
vives after extrapolation to the continuum limit. 

The difference in y^hirai jg ^^gg significant than it 
appears. All three systematic effects in HotQCD cal- 
culations discussed above as well as the change in 
WBC estimates on using the new experimental value 
of /if are again applicable. Even though the WBC 
find that the five quantities (rg , the new experimental 
value of /x, Mo, Mk' and M^) used to set the lattice 
scale give consistent estimates in the continuum limit, 
there can be systematic differences when working at a 
given a as is the case for the HotQCD data. Assuming 
these four ^ 5 MeV shifts add, as data indicate, the 
HotQCD estimate would be lowered to ^ 170 MeV in 
the continuum. This, however, still leaves a difference 
of about 24 MeV that needs to be resolved. 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



The good news is that the HotQCD simulations un- 
derway will allow us to extrapolate to the continuum 
limit and to the physical light quark mass and thus 
shed light on these issues. We hope to have these 
results within the next year. 
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